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I analyze an atomic Fermi gas with a planar p-wave interaction, motivated by the experimentally 
observed anisotropy in p-wave Feshbach resonances. An axial superfluid state is verified. A domain 
wall object is discovered to be a new topological defect of this superfluid and an explicit solution has 
been found. Gapless quasiparticles appear as bound states on the wall, dispersing in the continuum 
of reduced dimensions. Surprisingly, they are chiral, deeply related to fermion zero modes and 
anomalies in quantum chromodynamics. The chirality of the superfluid is manifested by a persistent 
anomalous mass current of atoms in the groundstate. This spectacular quantum phenomenon is a 
prediction for future experiments. 
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The p-wave Feshbach resonance has recently been 
demonstrated accessible and robust in atomic gases of 
6 Li and 40 K 0,0. Pioneer theoretical studies [Hi HQ 
point to interesting new properties associated with the 
resonance. For a simple model, one may consider a sys- 
tem of spinless fermionic atoms interacting through a p- 
wave potential, isotropic in the 3D orbital space. Such an 
interaction would enjoy the usual SO (3) spherical sym- 
metry in orbital space. However, the p-wave Feshbach 
resonances are split between the perpendicular and par- 
allel orbitals with respect to the magnetic field axis, due 
to the (anisotropic) magnetic dipole-dipole interaction. 
This is known both experimentally and theoretically 0. 
The energy splitting, in units of Bohr magneton, is signif- 
icantly greater than both the free atom Fermi energy [for 
a typical gas density] and the resonance width [which 
is narrow]. Therefore, it is perhaps best to distinguish 
the p XiV orbitals from the p 2 orbital and model the two 
p-wave resonances separately. 

The planar p-wave model Focus on the resonance in- 
volving the p XtV orbitals. I postulate a planar p-wave 
model for the atomic Fermi gas with the interaction 
parametrized by a single coupling constant g. Using a k 
and at as creation and annihilation operators for atoms 
with momentum k, the Hamiltonian reads 
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where Ck = j— ~ M is the energy spectrum of the free 
atom, measured from the chemical potential p, and 'f is 
the space volume {% = 1 unless explicitly restored). A 
convention is adopted: the boldfaced vector labels three 
components k = (k x , k y , k z ) and the arrow vector collects 
the x,y components only, k = (k x , k v ). When expanded 
in terms of the spherical harmonics, the interaction term 
contains only p x and p y orbitals as indicated by the arrow 
vector; it explicitly breaks the 3D rotational invariance 
down to a lower symmetry of SO (2) that rotates about 
the z-axis. Only L z of the three angular momenta is con- 
served. To compare with a physical system, the coupling 
constant g is related to the p-wave scattering length a\ 
and effective range ro by g ~ 2ivh ' J 2lT ° , with a\ < Q . 



As I shall show below, this atomic gas model predicts 
remarkable properties, including a new domain wall soli- 
ton and an anomalous quantum mass flow of atoms per- 
sisting in one direction along the wall. 

Axial superfluid For g < 0, the interaction term in 
Hamiltonian can be decoupled by introducing a pla- 
nar vector field $ via the Hubbard-Stratonovich trans- 
formation, standard in the path integral formalism of 
many-body theory. The field operator <£> q is identified 
by the quantum equation of motion with a composite, 
p-wave atom pair, 
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$ is the superfluid order parameter. This order parame- 
ter is a complex, two-component vector in orbital space. 
We parametrize it with four independent real variables 
as follows, 
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where 1 is a 2 x 2 identity matrix and 02 the second Pauli 
matrix. In literature, one encounters other popular pa- 
rameterizations, such as writing the order parameter as 
a rank-2 tensor or a sum of real and imaginary vectors. 
Whatever other representations are can be mapped, in 
a one-to-one correspondence manner, to the above form. 
An advantage of our parametrization is that its symme- 
try transformation property is made transparent. The 
matrices 1 and 02 are the generators of the U(l) phase 
and SO (2) orbital rotation transformations, respectively. 
So, i? represents the phase of the atom-pair wavefunction 
and if is the azimuthal angle (about the z-axis) of the 
order parameter as a planar vector in orbital space. For 
any finite modulus p, both symmetries are spontaneously 
broken down. The third angle variable, transforms 
X — * —X under the time reversal symmetry. 

For a homogeneous state of (spatially) uniform i>, the 
effective potential (or free energy if extended to finite 
temperature) density can be exactly calculated. I found 
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Notice that the potential does not depend on the phase 
$ and the orbital angle tp because of the symmetries. 
(Set t? = if = hereafter). Given a constant g, the ef- 
fective potential is minimized at a finite value p = po 
and at x = ±f,7r ± f (see Fig. [T] and Eq. ©). p is 
not universal, depending on the Fermi and ultra-violet 
cutoff momenta. The four minima of x are degenerate. 
A 7r phase identification, due to the U(l) phase sym- 
metry of the order parameter, reduces the four into two 
identical groups, each of two minima. So the angle vari- 
able x should be limited to half the scope of 2ir, say 
in [—f,^], to avoid the parametrization redundancy. I 
emphasize that the two minima x — + f an d X = — f 
are discrete due to time reversal symmetry. This will be 
of important topological implications later. Then, the 
(mean-field) states are the (axial) p x ± ip y superfluid: 
$0 = "y|(±i) f° r X — if j respectively. One of these two 
states will become the groundstate, breaking the sym- 
metry spontaneously. My finding agrees to the results of 
other different p-wave models |j, |a] ■ The state carries a 




FIG. 1: Effective potential V B a of Eq. I^J per unit volume, 
with coordinates [<j> x , i(j> y ] — [pcosx, — psinx] in order param- 
eter space and (?/(2-7r) 3 = 0.5. The energy and length units 
are the Fermi energy (equal to p) and the Fermi wavelength 



k F , respectively. 



total angular momentum, estimated to be L* otal 
with TVn the number of atoms in the condensed pairs. Nq 
is directly accessible in current experiments. 

Domain Walls The degeneracy of the x = if states 
motivates us to look for a spatially nonuniform, soli- 
ton solution that intervenes the two minima and varies 
smoothly from one to another in real space. Such a con- 
figuration is topologically stable since the two minima are 
related by the discrete time reversal Zi symmetry, not by 
any continuous symmetries such as orbital rotation and 
phase transformation. 

To uncover the soliton granted by the topology of bro- 
ken time-reversal symmetry, we consider temporally sta- 



tionary, spatially nonuniform configurations of x( r ) with 
r the space coordinate while keeping the modulus con- 
stant. The order parameter then becomes 
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where l>(r) = ^ q eiq r *q- 

Our goal is to find an effective theory for the angle 
field x( r ) that retains both spatial derivative terms and 
potential terms. Such an effective theory is obtained per- 
turbatively through the Wilsonian renormalization ap- 
proach, by integrating out fermionic atoms at high ener- 
gies (measured from the Fermi level) while keeping low 
energy fermions within a thin shell of thickness, say k, 
along the Fermi sphere in the momentum space. 1 have 
found the effective free energy functional of x( r ); 
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The characteristic length £, emergent in the low energy 
effective theory, is the coherence scale on which the or- 
der parameter varies. A perturbative calculation deter- 
mines £ = -^4- x Cf(ir-, tM where An is the maxi- 
mum energy gap of quasiparticle excitations. C$ and 
are dimensionless positive constants depending on the 
infrared and ultraviolet cutoffs (k and A), and kF is the 
Fermi wavevector hp /(2m) = p. [Detailed derivations, 
including coefficients C$ and C{, will be given elsewhere.] 
From either this perturbative calculation or the general 
Ginzburg-Landau phenomenology, the order parameter 
should vary slowly compared with the Fermi wavelength, 
i.e., l/(£,kp) <C f. The effective theory is invariant under 
a Z2 time-reversal symmetry, x —* ~X- 

The effective theory happens to be a variant of the 
well-known sine-Gordon Lagrangian in the 3D Euclidean 
space. It is known to possess (nonuniform) soliton solu- 
tions 0]. Applying the textbook method, we find soliton 
and anti-soliton configurations (say, constant in y, z di- 
rections) , 
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X±(r) = iarctanl tanh(a;/£) 



where '±' are the two topological charges (or winding 
numbers). For either configuration, x( r ) changes sign 
when x varies from —00 to +00, developing a wall that 
separates two domains of opposite angular momentum. 
The domain wall is centered at x = 0, a position however 
set arbitrarily for simplicity. The thickness of the wall is 
characterized by £. Its energy per unit area is calculated 
to be Po/(^*0; proportional to A /^ 2 (recall An the 
maximum energy gap). 

Gapless chiral fermion bound states What happens 
to the quasiparticle states when the order parameter 
develops a soliton defect? Consider an order parame- 
ter configuration in real space depicted in Fig. x 
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varies in a characteristic distance of £, far greater than 
the Fermi wavelength kp 1 . The fermionic atoms view 
the soliton defect as a fairly flat, "classical" off-diagonal 
potential that hybridizes the particle and hole states of 
atoms. So we can loosely speak of three different re- 
gions along the x-axis — left (p x — ip y ), domain wall (p x 
state), and right (p x + ip y ) boxes. Now imagine a mo- 
mentum space for each 'macroscopic box', each beginning 
with free fermionic atoms labeled by 3-component defi- 
nite momentum k. An important observation is that the 
off-diagonal pairing potential seen by the atomic states 
in the fcj,-axis direction must be vanishingly small in the 
domain wall box (exact zero at the center of the box) 
but maximizes at infinities in the left and right boxes. 
Therefore, like the Caroli-de Gennes-Matricon state 0] 
in the vortex core of a superconductor where the order 
parameter vanishes, we expect to see zero energy states 
peaked at the domain wall and bound by the rising p y - 
wave component of the gap potential at far away from 
the wall. Unlike the vortex core states, the domain wall 
fcrmions are only bound in the direction perpendicular 
to the wall (x-direction) and disperses like in continuum 
in the parallel directions (Fig. 2). 
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FIG. 2: Heuristic illustration of gapless chiral quasiparticlc 
bound states: (a) momentum-space energy gap in the axial 
Px+ipy superfluid state; (b) gap in thep x state; (c) the soliton 
defect in the real space with the angular momentum L — 
— i$* x $ = e z po sm Xi (d) the gapless bound state profile 
with \<& y \ as the "off-diagonal" potential ('■••')• The arrow 
in (c) indicates the direction of anomalous mass flow. 



where k|| = (0, k y ,k z ) is a momentum parallel to the 
domain wall, and n a quantum number from the quanti- 
zation of x-direction. The unitarity of the transformation 
requires that / c2 3 r[u„ k|| (r)«;, k , (r) + u nk|| (r)u*, k , (r)] = 

^nn' ^k||,kj| ■ The unbroken symmetries of translation par- 
allel to the domain wall simplify the eigenstates, 

/ \ / \ e ik ii' r 

(u nk|| (r) , v nk[l (r) J = [u nk[[ (x) , u„ k|| fx) J — — . 

We are primarily interested in finding the low energy 
quasiparticle states that have a parallel momentum in or- 
der of |ky | ~ kF (Fermi wavevector) and vary slowly in x 
direction in real space. For £kp ^> 1, the eigenstates are 
dictated by the following Bogoliubov-de Gennes equa- 
tions 

-ip (cosxd x - kySmxjVn^ = (-E„ k|| - e k|| )u nk|| , 
-ip (cosxd x + fcj,smx)u nk|| = (-EWk y + ek y )«nk|| , 

where x = X+( r ) [Eq. 0] and po is a constant related 
to the maximum gap by pokp = Ao. 

The above equations have symmetries. While revers- 
ing the sign of both x an( A k v simultaneously or revers- 
ing the sign of k z independently, the eigenvalues are in- 
variant. One can also prove that given any eigenstate 
(w„kii ! v n kii ) with eigenvalue S nk|| , there always exists an- 
other eigenstate, (— u nk , u nk ), with an eigenvalue of op- 
posite sign, — E n k , where k|| = (0, — k y , k z ). Such a dou- 
bling of eigenstates are not surprising in a superfluid sys- 
tem that had originally enjoyed time-reversal and parity 
symmetries before undergoing superfluid. However, as I 
shall show later, the phenomenon of doubling eigenstates 
breaks down when chiral fermions appear. There, the 
physical eigenstates are no longer symmetrically available 
upon reversing the sign of momentum, say k y — > — k y . 

The reduced ID eigenvalue problem resembles, but dif- 
fers in fundamental ways from, the Dirac equation of one 
spatial dimension analyzed by Jackiew and Rebbi \Ti\ 
who found the chiral zero-energy bound state (zero 
fermion mode). In our case, we would expect that their 
zero-energy mode translates into gapless chiral quasipar- 
ticle bound states (chiral fermions). Indeed, focusing on 
the gapless branch (n — 0) , I found two branches of such 
states. They are 
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Such heuristic argument can be made rigorous by an- 
alyzing the Hamiltonian QJ in the axial state with the 
domain wall defect. The (mean- field) Hamiltonian can 
be diagonalized by using the Bogoliubov transformation 
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with '~' meaning a proper normalization to be achieved 
yet. Clearly, the E + and E~ branches are only normal- 
izable (so become physical bound states) for k y > and 
k y < 0, respectively. The superscript '±' is best inter- 
preted as a sign of positive or negative chirality for the 
existing gapless fermion modes. For either chirality the 
energy S nk| . can be positive or negative, depending on 
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the value of = kjj/ (2m)— fx. The doubling of eigenval- 
ues is removed by the presence of chirality, with only one 
mode being physical bound state, supporting the claim 
made early. For gapped branches (denoted with n > 0), 
I would speculate, based on the study of domain wall 
quasiparticle excitations in 3 He 0, that bound states 
exist for both positive and negative momentum k y (hence 
not chiral) and the doubling should restore. 

Anomalous quantum mass flow With the eigenfunc- 
tions diagonalizing the superfluid Hamiltonian, we can 
directly calculate the atom occupation number in the 
state vector space, specified by the quantum numbers 
nky). I found that the atom occupation number per 
state is N n -^.. = j dx\v n ^(x)\ 2 , derived from the to- 
tal atom number N — f d 3 r(a^ (r)a(r)) = ^ nk A^ky ■ 
Clearly, the discrete quantum number n takes the role 
of k x for a homogeneous superfluid state (where the do- 
main wall defect is absent). Finite temperature effects, 
which are not included in this paper, will modify the 
above result through extra terms and factors dependent 
on the Fermi-Dirac distribution function. The occupa- 
tion weight is essentially given by the eigenfunction v 
alone, consistent with the result established for the homo- 
geneous superconductivity. The physical normalizability 
constrains that non-zero Vquu exist only for negative chi- 
rality of k y < (Eq. @). We then must conclude that 
the atoms of quantum number n = only occupy half the 
reduced 2D k y -k z momentum plane (parallel to the do- 
main wall), with a spectral weight |£>ok||| 2 - The another 
half space is exactly empty for the gapless states in the 
superfluid groundstate (zero temperature). Bear in mind 
that the n > states are different. This is the origin of 
an anomalous chiral mass flow of atoms. 

A direct calculation of the mass current of atoms veri- 
fies the above intuitive speculation. I found that the chi- 
ral gapless bound states give rise to a spectacular anoma- 
lous current of atoms equal to 

hk 3 

j = —-r-7re v , (H restored for clarity) 

OTT z 

per unit length in z-axis (perpendicular to the flow). The 
current persists without an additional external field. Its 
direction is set by the topological charge of the soliton 
defect. For an anti-soliton, the anomalous mass flow re- 
verses direction. In real experiments that the magnetic 
field sets the z-axis (H || e z ), the prediction is that the 
anomalous mass current flows parallel to the domain wall 
in the direction of VL Z x e z where VL Z is taken as the gra- 
dient of the angular momentum L z at the wall (Fig. 2c). 
The mass flow returns at the boundaries of the atomic 
trap which constitute an anti-soliton (anti-domain wall) . 
The total mass current, therefore, does not violate con- 
servation law. 

At finite temperature, excitations with positive energy 
will proliferate those unoccupied states in the empty half 
space and/or deplete the occupied states. The finite 



temperature effect is then to create a counter mass flow 
in opposite direction. I conjecture that the anomalous 
mass flow should decrease with increasing temperature 
and should be strongest in magnitude for temperatures 
(times the Boltzmann constant) below the level splitting 
of quasiparticle bound states, characteristically given by 

Discussion The anomalous current is reminiscent of 
that in the 3 He A-phase [HE1E3 and in the PbTe semi- 
conductor |lfil |. but differs in nature. First, the domain 
wall defect I discovered is a new kind different than the 
"twist" texture of Ref. [HHQ for 3 He-A. In the lat- 
ter, the angular momentum L sweeps its direction in real 
space with a unit magnitude. This domain wall is how- 
ever that L is fixed in the z-axis direction but changes 
sign and magnitude across the wall (Fig. |2t). Second, 
treating the pairing potential scmi-classically, the quasi- 
particles see a gapless line (in the k y -k z plane of Fig. EJ^) 
in the domain wall region in this case as opposed to nodal 
points everywhere in the 3 He-A case. These features, 
among others, are new for the atomic gas of anisot ropic p- 
wave Feshbach resonances. The 3 He-A results [l2lll3.ll4T ] 
do not directly apply. 

This domain wall defect seems interesting and new for 
the experiment of atomic Fermi gases, currently develop- 
ing rapidly. I speculate that superfluid domains of oppo- 
site angular momentum may be easier to realize than a 
homogeneous (axial) superfluid, since the latter requires 
rotating the whole gas (or other means) for a net in- 
put of macroscopic angular momentum. I am not aware 
whether the anomalous current predicted for the 3 He-A 
has been observed experimentally. It would be first if 
the p-wave resonant atomic gas shows this spectacular 
quantum anomaly. That would be also of interest to the 
study of lattice quantum chromodynamics attempting to 
simulate chiral quark fields by fermion zero-energy mode 
bound at a domain wall llal. 
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